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Abstract 

The multifractal detrended fluctuation analysis technique is employed to ana¬ 
lyze the time series of gold consumer price index (CPI) and the market trend 
of three world’s highest gold consuming countries, namely China, India and 
Turkey for the period: 1993July 2013. Various multifractal variables, such as 
the generalized Hurst exponent, the multifractal exponent and the singularity 
spectrum, are calculated and the results are fitted to the generalized binomial 
multifractal (CBM) series that consists of only two parameters. Special empha¬ 
sis is given to identify the possible source(s) of multifractality in these series. 

Our analysis shows that the CPI series and all three market series are of multi¬ 
fractal nature. The origin of multifractality for the CPI time series and Indian 
market series is found due to a long-range time correlation, whereas it is mostly 
due to the fat-tailed probability distributions of the values for the Chinese and 
Turkey markets. The CBM model series more or less describes all the time 
series analyzed here. 
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Keywords: Multifractality; Detrended fluctuation analysis; Consumer price index; 

Long-range time correlation; Generalized binomial multifractal model 


1 Introduction 

In general, a fractal is a rough or fragmented geometrical shape that can be subdi¬ 
vided into parts, each of which is (at least approximately) a reduced-size copy of the 
whole. A fractal system is usually described by a scale invariant parameter called 
fractal dimension [T]. Many fractals arising in nature have a far more complex scaling 
relation than simple fractals and require a set of parameters to specify such objects 
that are known as multifractals. Several approaches have so far been developed and 
applied to explore the of fractal properties. For instance, the rescaled adjusted range 
analysis method was introduced by Hurst Eli (see also 0), which he himself applied 
the method to his hydrological study. Due to the difficulty of the rescaled analysis 
in capturing long-range correlations of nonstationary series, Peng et al. [5] proposed 
an alternative approach to analyze the DNA sequences which is known as detrended 
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fluctuation analysis (DFA). Although the DFA method is widely used to determine 
monofractal scaling properties, it cannot properly describe multi-scale and fractal sub¬ 
sets of time series data. One of the simplest type of multifractal analysis has been 
developed based on the standard partition function multifractal formalism (HE]. This 
is a highly successful method for the multifractal characterization of normalized and 
stationary measures, but it does not give the correct result for nonstationary time se¬ 
ries. Based on a generalization of the DFA method, Kantelhardt et al. [7] introduced 
the multifractal detrended fluctuation analysis (MF-DFA) for the multifractal charac¬ 
terization of nonstationary time series. As a remarkable powerful technique, MF-DFA 
has so far been applied to various fields of stochastic analysis, for instance, in mar¬ 
kets return analysis P P UHl EH E21 E31 EH E5], in geophysics [HI EH [201 ESI EO], in 
biophysics [21112211231121], and also in various branches of basics and applied physics 
[OslESlETlEElEnilao]. 

The study of financial time series has been the focus of intense research by the 
physics community in the last several years. Nowadays, there are some excellent 
compilations available on this subject, e.g. [SI11321133]) 1’^®^ ^o cite some of them. The 
prime objective of this kind of analysis is to characterize the statistical properties of 
the time series with the hope that a better understanding of the underlying dynamics 
could provide useful information to create new models. Moreover, such knowledge 
might be crucial to tackle relevant problems in finance, such as risk management or 
the design of optimal portfolios. Henceforth our discussion will be restricted to the 
time series analysis of gold market in China, Indian, Turkey and the global consumer 
price index (CPI)0. 

Gold being one of the most precious metals is always considered as the safest 
investment. Presently the fluctuations of gold market seem quite confusing even to 
the regular traders, and it becomes almost impossible to predict its accurate rise or 
fall. As we know, over the last 2/3 years gold price increases so rapidly that the 
gold price nowadays is approximately double of the average rate of 2010. The market 
gains its highest value of about 1900 USD/ounce in the year 2011, and the recent 
value is about 1300 USD/ounce. Moreover, the day-to-day variation of the market is 
also quite remarkable during the last few years. It is now believed that gold is not 
a commodity anymore, rather a currency which always maintain an inverse relation 
with the US economy. So it is quite possible that a part of the change in gold price 
is really just a reflection of a change in the value of the US dollar. Sometimes such 
change is insignificant and often the opposite is true. Whatever may be the reason, 
the dynamical nature of the gold market is quite complex, and one needs to study 
the time series of gold price from all possible directions in order to understand the 
underlying mechanism. 

In this article we apply the MF-DFA technique to characterize the time series of 
the gold CPI and the gold market in China, India and Turkey during the period 1993- 
July 2013. According to the World Gold Council, these three countries are the world’s 
major gold consuming countries with a combined consumption of about 70% of the 
total demand. The individual consumptions of these countries are: China 33%, India 
28% and Turkey 9%. Hence, it is expected that the CPI of gold is mainly governed by 

consumer price index (CPI) is an estimate as to the price level of consumer goods and services 
in an economy which is used as a way to estimate changes in prices and inflation. A CPI takes a 
certain basket of common goods and services, for instance a gallon of gasoline and diesel fuel or an 
ounce of gold, and tracks the changes in the prices of that basket of goods over time. The gold 
CPI, according to the World Gold Council [33], is composed of the 5 largest gold consuming country 
currencies, ranked by and weighted by 3 year average gold demand. 
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the markets of these three countries. In order to visualize the recent market pattern, 
we separately analyze the series of about the last three years period-from 2010 to 
June 2013. Various parameters related to multifractality of the studied time series 
are computed and are htted to the two-parameter generalized binomial multifractal 
(GBM) model [H [12]. Special emphasis is given to identify the probable origin(s) 
of mult ifr act ality in these series. For this purpose we analyze a randomly shuffled 
series and a surrogate series corresponding to each of the original series. The article is 
organized as follows: in Section 2 the MF-DFA methodology is described along with 
the outlines of the GBM model. In Section 3 we describe the data and the results of 
our analysis, and the article is summarized in Section 4. 

2 MF-DFA methodology 

Though nowadays the MF-DFA technique has become a standard tool of time series 
analysis, for the sake of completeness we provide in this section a brief description of 
the method which is followed by the outlines of the GBM model used to compare the 
empirical data. 

Let {xk '■ A; = 1, 2,..., V} be a time series of length N. The MF-DFA procedure 
consists of the following hve steps: 

Step 1: Determine the profile 


^(*) = * = 1, 2, ...,V, (1) 

k=l 

where (x) = (l/N) Yl!k=i^k is the mean value of the analyzed time series. 

Step 2: Divide the profile Y[i) into Ng = int{N/s) non-overlapping segments of equal 
length s. One has to choose the s value depending upon the length of the series. 
In the case, the length N is not a multiple of the considered time scale s, the same 
dividing procedure is repeated starting from the opposite end of the series. Hence, in 
order not to disregard any part of the series, usually altogether 2Ns segments of equal 
length are obtained. 

Step 3: Galculate the local trend for each of the 2Ns segments. This is done by a 
least-square fit of the segments (or subseries). Linear, quadratic, cubic or even higher 
order polynomial may be used to detrend the series, and accordingly the procedure 
is said to be the MF-DFAl, MF-DFA2, MF-DFA3, ... analysis. Let i/p be the best 
htted polynomial to an arbitrary segment p of the series. Then determine the variance 

F\p, s) = \Y. - 1)* + 'I - (2) 

i=l 

for p = 1,... Ng^ and for p = -|- 1,..., 2Ng it is given as. 


F\p, s) = -Y,{y[N-{p- Ng)s + z] - yp{{)y 


2=1 


Step Dehne the gth order MF-DFA huctuation function 


2Ns 


pys) = 


2N. 




p=i 


j' 


(3) 


(4) 
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for all g 7^ 0 and for g = 0 the above definition is modified to the following form 


Fg{s) = exp 


1 

4 ^ 


^ln[F2(p,s)] 


p=i 


( 5 ) 


Step 5: Then the scaling behavior of the fluctuation functions is examined for several 
different values of the exponent g. 

If the series {xfc} possess long-range (power-law) correlation, Fq{s) for large values 
of s would follow a power-law type of scaling relation like 

F,{s) ~ s'-!'*. (6) 


In general, the exponent h{q) depends on g and is known as the generalized Hurst 
exponent. For a stationary time series h{2) = Ff - the well known Hurst exponent 
|16j . For a monofractal series on the other hand, h{q) is independent of g, since the 
variance F‘^{p,s) is identical for all the subseries and hence Eqns. (jl]) and (jH]) yield 
identical values for all g. Note that the fluctuation function Fq{s) can be dehned only 
for s > m -|- 2, where m is the order of the detrending polynomial. Moreover, Fq{s) 
is statistically unstable for very large s (> N/4). If small and large fluctuations scale 
differently, there will be a significant dependence of h{q) on g. For positive values 
of g, Fq{s) will be dominated by the large variance which corresponds to the large 
deviations from the detrending polynomial, whereas for negative values of g, major 
contributions in Fq{s) arise form small fluctuations from the detrending polynomial. 
Thus, for positive/negative values of g, h{q) describes the scaling behavior of the 
segment with large/small fluctuations. 

One can easily relate the h{q) exponent with the standard multifractal exponent, 
such as the multifractal (mass) exponent T(g). Consider that the series {xk} is a 
stationary and normalized one. Then the detrending procedure in step 3 of the MF- 
DFA methodology is not required, and the variance of such series is given by 

FUp,s) = {Y{ps)-Y[{p-l)s]y. (7) 


Then the fluctuation function and its scaling law are given by 

( ^ 2Ns '\ 

Now if we assume that the length of the series N is an integer multiple of the scale s, 
then the above relation can be rewritten as, 

N/s 

^|K(ps)-y|(p-l)s]r~ii’'*l«l-‘. (9) 

P=1 


In the above relation the term under | • | is nothing but sum a of {xk\ within an 
arbitrary pth segment of length s. In the standard theory of multifractals it is known 
as the box probability V{s,p) for the series Xk- Hence, 

ps 

F{p,s)= ^ Xk = Y{ps) -Y{{p-l)s). (10) 

/c=(p—l)s+l 
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The multifractal scaling exponent T{q) is defined via the partition function Zp{s) 

N/s 

zp(i,)s^iT’(p,i,)r~ V’>, (11) 

p=i 

where g is a real parameter. From Eqns. ©-(HI]) it is clear that the multifractal 
exponent r(g) is related to h{q) through the following relation: 

r(g) = qh{q) - 1. (12) 

Knowing r(g) one can calculate the most important parameter of a multifractal anal¬ 
ysis - the multifractal singularity spectrum (also called the spectral function) /(a) 
which is related to r(g) through a Legendre transformation [H |35]: a = dT{q)/dq, 
and is given by 

f{a) = qa-T{q). (13) 

Here a is the singularity strength or Holder exponent. Importance of the singularity 
spectrum in the theory of multifractals is that, the width of the spectrum is a direct 
measure of the degrees of multifractality. For a monofractal object it turns out to be 
a delta function at the corresponding a. 


The binomial multifractal series of numbers is defined as, 


with A: = 1,..., iV, n{k) is the number of digits equal to 1 in the binary representation 
of the index k and 0.5 < a < 1 is a parameter. A generalization of this series with 
two positive valued parameters (say, a and b = (1 — a)) is given as. 


_ ^n(k—l)^nma,x—n{k—l) 


(15) 


The above generalized form of the binomial multifractal series (ITT)) has been used in 
[TB] to characterize the river runoff data. In the text it is said to be the generalized 
binomial multifractal (GBM) series. With this generalization one can easily derive 
the expressions for the multifractal observables h{q) and T{q) of the following form: 


h{q) 
and r(g) 


1 ln[a'^ -|- 6'^] 
q g In 2 
+ h^] 
K2 ’ 


(16) 

(17) 


Note that for a = b, h{q) = — In a/In 2 is independent of q and r(g) ~ gin a/In 2 is a 
linear function of g, i.e. for a = b the series flT5|l reduces to a monofractal series. Then 
obviously the quantity |a — 6| is a measure of the strength of multifractality. Quan¬ 
titatively, the strength parameter Ack is defined by the difference of the asymptotic 
values of h(g), i.e., Ack = h{—oo) — h{+oo) = (Ina — In6)/In2, which is nothing but 
the width of the singularity spectrum at /(ck) = 0. 


3 Results and discussion 

The gold market data used in this article are taken from the database of the World 
Gold Gouncil [3l]. The logarithmic difference between two successive trading days. 
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Figure 1: (a) Time series of the gold consumer price index and the corresponding 
return for the time period 1993-July 2013, (b) same as (a) but for the time period 
2010-July 2013. 


also known as return: R{T) = In P(T + 1) — InP(T), is used as the original series for 
this analysis. To illustrate the nature of the series we show in Fig. [U^a) the CPI time 
series (upper panel) along with its returns (lower panel). The last section of these 
series shown in blue color are for the period of 2010-July 2013 and is magnihed in 
Fig. [11(b). In the text the hrst section, i.e. from 1993 to 2010, of the series is identihed 
as series of period-I and the last section (from 2010 to July 2013) is identihed as period- 
II. Note that, the time periods are explicitly mentioned in the diagrams. Figure [1] 
demonstrates how the global gold market evolves with time. In Fig. [2] we show the 
log-log plots of the cumulative distribution function (CDF) of normalized returns of 
gold price indices for (a) the CPI and (b) the market series in China. The CDFs 
for the other series analyzed here are more or less similar to [2](b) and hence are not 
shown here. However, the tail exponent, {(), gives the critical order of divergence of 
moments, is extracted from power-law regression for all the series. The ( values are 
all found to vary from 3.04 to 3.5 for the long series, whereas for the series of period-II 
(2010-2013) the values are as high as 5.0 (for the CPI series). 

The hrst indication for a power-law type of CDF of market return huctuations can 
be traced back to [3S] and explicitly in stock markets [5H1 [3Z| the power-law is found 
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Figure 2 : Cumulative distribution function of the normalized returns of gold price 
indices for three different time periods: (a) global consnmer price index and (b) time 
series of China. 

to be inverse cubic. Gabaix et ah [39] have given a theoretical interpretation of the 
above empirical observations. In the case of gold market returns for the period 1968 
- 2010 the tail exponent is estimated to be ~ 3, i.e. the CDF follows inverse cnbic 
power-law [ID] . The effect of faster departure from the inverse cubic power-law in the 
more recent hnancial data, similar to the one shown in the present manuscript for the 
gold price time series in the period 2010-2013, is also shown in [TTlIT3] . Recently, 
Rak et ah mi indicate a possible more general applicability of the concept of Gabaix 
et al. [39] to the situations when the price flnctnations depart from the inverse cubic 
power-law. 

Long term correlated records {xi : i = 1,2,... N} with zero mean and nnit vari¬ 
ance are characterized by an antocorrelation fnnction Cx{s) = {xiXi+s) = — 

s) If cire nncorrelated, C'(s) is zero for s > 0. Short-range correla¬ 

tions of {xi\ are described by C'(s) declining exponentially: G(s) ~ exp(s/sx) with a 
decay time Sx. For long-range correlations G(s) declines as a power-law: G(s) ~ , 

where the correlation exponent 7 is between 0 and 1. A direct calculation of C{s) is 
usnally not preferable becanse of the nnderlying trends of the series {xi} and its noise 
factor. To test for long-term correlations, it is nsefnl to employ the detrended flnctn- 
ation analysis (DFA) in order to extract the Hurst [H) exponent, which is related to 
the 7 exponent as: H = 1 — 7/2 [D]. However, as a preliminary estimate of the natnre 
of correlation present in the gold price time series data analyzed here, we evaluate the 
autocorrelation fnnctions for the series. In Fig. [3] we illustrate the autocorrelation 
fnnctions for the consnmer price index time series as well as for the Indian market 
time series. We observe that the daily returns of gold indices, irrespective of the 
market, are nncorrelated (Fig. [3](a)). Bnt the maxima and minima of the individnal 
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Figure 3: Autocorrelation functions for the gold consumer price index (CPI) series 
compared with the price time series in Indian market, (a) Includes total series for 
the period 1993 - July 2013, which is totally uncorrelated, (b)-(f) are drawn for the 
maxima/minima returns in windows of length 5 for several different time periods as 
mentioned in the diagrams. Notice the signature of long-range correlation in (b)-(d). 
The correlation function in (b) for a period of about 6 years captures a very weak 
correlation, whereas (e) and (f) do not possesses any noticeable correlation pattern. 


market returns over a reasonable time duration (at least 5-6 years) in windows of 
length R = 5 days (1 trading week) separately exhibit long-range correlations (Fig. 
[3](b)-(f)). Therefore, we do not expect long-range correlation in the series of period 
2010-2013. Note that the series of maximum (minimum) returns corresponding to 
each of the original ones are constructed by selecting maximal (minimal) values of 
daily returns over each interval of length R. Obviously, both sub-sequences contain R 
times fewer points than the original ones. Another importance of the market return 
autocorrelation function is that, it can serve as a proxy for market liquidity, though 
the later is an elusive concept, and hence is not a direct observable. Moreover, market 
liquidity has several dimensions of measurement (e.g., see BSI). So a rigorous estimate 
of the gold market liquidity is beyond the scope of the present manuscript, which will 
otherwise be conhned to the short/long-range correlations of gold market returns and 
its (multi)fractal nature. 

We calculate the MF-DFA2 fluctuation functions Fq{s) for all the series under 
consideration and over a wide range of the exponent q : -10 to -|-10 in steps of 0.5. 

















































The scaling behaviors of some of the fluctuation functions calculated from the CPI 
series are illustrated in Fig. HI Separate diagrams are drawn for the total series, 
period-I and period-II. Here the time scale s is varied from 6 to iV/5, where N is the 
series length. In the same hgure we also show the Fq{s) functions generated from the 
randomly shuffled series (middle panel) as well as from the surrogated series (lower 
panel) corresponding to the original ones. The fluctuation functions calculated for 
the other series under consideration also follow almost identical scale dependence as 
shown in the figure, and hence those diagrams are not included here. It is noticed that 
the multifractal results reasonably vary between the first and second order detrending, 
but the change is almost insignificant between the second and third or higher order 
detrending. Hence, the quadratic (second order) detrending is considered for this 
analysis. It is clear from Fig. [T]that the scaling behavior of Fq{s) for all three time 
intervals as well as their shuffled and surrogated counterparts nicely follow the scaling- 
law ([6|) for s > 10. 

It is a known fact that there are two different types of multifractality may exist 
in a time series data, namely (i) mult ifr act ality due to long-range time correlations of 
the small and large fluctuations and (ii) multifractality due to a fat-tailed probability 
distribution function of the values in the series. The first kind of multifractality can be 
removed by random shuffling of the given series and the corresponding shuffled series 
will exhibit monofractal scaling. Obviously, the probability distribution will not alter 
by random shuffling and hence the multifractality of the second kind will remain 
intact. If a given series contains both kinds of multifractality, the corresponding 
shuffled series will exhibit weaker multifractality than the original one. On the other 
hand, the surrogate (phase randomization) analysis is an empirical technique of testing 
nonlinearity for a time series. The aim is to test whether the dynamics are consistent 
with linearly hltered noise or a nonlinear dynamical system gSliZ]. The basic idea of 
the surrogate data method is to first specify some kind of linear stochastic process that 
mimics “linear properties” of the original data. If the predictions (statistics) of the 
original data are significantly different from those of surrogate series, we may consider 
the presence of some higher order temporal correlations, i.e., the presence of dynamic 
nonlinearities. In this analysis we use the the null hypothesis, amplitude-adjusted 
Fourier transform (AAFT) algorithm [IS] to generate the surrogate series. 

The generalized Hurst exponent h{q) is extracted from straight line £t to the log-log 
data of Fq{s) versus s. We fit straight line in the 50 < s < 800 region for the total and 
the period-I series, whereas the range chosen for the period-H series is 10 < s < 60. 
Note that, the fluctuation functions are found to obey the scaling law (|HD better in 
the mentioned scale regions. The order (q) dependence of the h{q) exponents is shown 
in Fig. [5] for the CPI series as well as for the price time series in China, India and 
Turkey. The lines in the diagrams represent the best fits of the GBM model (whenever 
possible) to the data points. The best fit parameter values are given in Table [T] 
The missing values in the table (and line in the figure) denote that the GBM model 
cannot describe the corresponding time series. The figure shows that, irrespective of 
the series time duration, all the h{q) spectra representing the original series are order 
dependent, though the dependency is weaker for the series of period-H. Moreover, the 
spectra obtained from the original series are more or less fitted to the GBM model, 
yielding Eqn. flTB]) . hOWEVER, in some cases there is noticeable deviation between 
the empirical values and the model around g ~ 0 regions. The shuffled series, on 
the other hand, produces h{q) spectra that are quite different from their original 
counterparts. In all cases but the GPI series of period-H, the total and period-I series 
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Table 1: The best fitted parameter (a and b) values of the generalized binomial mul¬ 
tifractal series flT5|l and the multifractal strength parameter: Aa = (Ina — In 6)/In2. 
The missing data indicate that the binomial series is unable to describe the corre- 
sponding time series. 


Origin 

Series period 

Series type 

a 

b 

Aa 



Original 

0.575±0.008 

0.904±0.013 

0.652T0.020 


1993 to July 2013 

Shuffled 

0.639±0.012 

0.786±0.015 

0.297±0.027 



Surrogated 

0.648±0.011 

0.802±0.014 

0.306±0.025 



Original 

0.576±0.008 

0.917T0.013 

0.670T0.021 

CPI 

1993 to 2010 

Shuffled 

0.645±0.011 

0.805±0.014 

0.319±0.025 



Surrogated 

0.651±0.011 

0.826±0.013 

0.344±0.023 



Original 

0.632±0.011 

0.805±0.014 

0.351T0.025 


2010 to July 2013 

Shuffled 

— 

— 

— 



Surrogated 

— 

— 

— 



Original 

0.572T0.008 

0.933T0.014 

0.706T0.021 


1993 to July 2013 

Shuffled 

— 

— 

— 



Surrogated 

— 

— 

— 



Original 

0.568±0.008 

0.938±0.013 

0.724±0.021 

China 

1993 to 2010 

Shuffled 

— 

— 

— 



Surrogated 

— 

— 

— 



Original 

0.640±0.011 

0.816T0.013 

0.351T0.023 


2010 to July 2013 

Shuffled 

0.639±0.009 

0.852±0.012 

0.416±0.021 



Surrogated 

— 

— 

— 



Original 

0.585±0.008 

0.929±0.013 

0.667±0.020 


1993 to July 2013 

Shuffled 

0.635±0.013 

0.678±0.016 

0.291±0.028 



Surrogated 

0.662±0.012 

0.798±0.014 

0.272T0.026 



Original 

0.593±0.009 

0.946T0.014 

0.674T0.021 

India 

1993 to 2010 

Shuffled 

0.654±0.013 

0.779±0.016 

0.252T0.029 



Surrogated 

0.659±0.013 

0.790±0.015 

0.262T0.027 



Original 

0.630±0.011 

0.816T0.013 

0.374±0.024 


2010 to July 2013 

Shuffled 

0.623±0.010 

0.814±0.013 

0.385±0.024 



Surrogated 

0.636±0.010 

0.824±0.013 

0.372T0.019 



Original 

0.528±0.010 

0.853T0.016 

0.629T0.026 


1993 to July 2013 

Shuffled 

— 

— 

— 



Surrogated 

— 

— 

— 



Original 

0.541±0.011 

0.862±0.017 

0.671±0.028 

Turkey 

1993 to 2010 

Shuffled 

— 

— 

— 



Surrogated 

— 

— 

— 



Original 

0.618±0.012 

0.771T0.015 

0.319T0.028 


2010 to July 2013 

Shuffled 

0.663±0.017 

0.754±0.020 

0.185±0.038 



Surrogated 

— 

— 

— 
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of China and Turkey, we observe a very weak order dependence of h{q) and the spectra 
are well described by the GBM model. In view of our previous discussion, the origin of 
multifractality in those series are both long-range temporal correlation and fat-tailed 
probability distribution of the variables. In the exceptional cases, the shuffled series 
produces h{q) spectra that are apparently consistent with an uncorrelated multifractal 
series with power-law distribution function: P{x) = for 1 < x < cxo and 

a > 0. For the above series the generalized Hurst exponent is given as 0 : K(l) ~ 1/9 
for g > a and ~ 1/a for q < a, whereas r(g) = q/a — 1 for g < a and r(g) = 0 
for g > a. Note that the multifractality of the above kind is originated purely from 
a fat-tailed probability distribution function. According to refs. [1^1 [50] . the CPI 
series for period-II, and the total and period-I series of China and Turkey, may be 
identihed as bifractals instead of multifractals. Our results on h{q) spectra apparently 
contradict the results of ref. M. where it was possible to remove the multifractality 
of the original series by random shuffling. In a similar analysis by Ghosh et al. [51] 
a g-dependent h(g) spectrum is observed for the shuffled series of gold indices. But 
one should keep in mind that the series time period in those analysis are different, 
and the recent rapid changing values of gold indices may moderately change the 
results. In our analysis the surrogate data generated by the AAFT procedure also 
result in h{q) spectra that show significant amount of multifractality. The conjectures 
will be more transparent when we will discuss the singularity spectrum. It is to 
be noted that, in order to minimize the statistical error in the multifractal variables 
corresponding to the shuffled/surrogated series, we take an average of h(g) values over 
10 shuffled/surrogated series corresponding to each of the original series. 

The multifractal r(g) exponents are plotted against g in Fig. [6] for all the series 
under consideration. The lines in this hgure also represent the binomial series, yielding 
Eqn. fITTD . with the same sets of parameters that are used optimiz the h{q) versus g 
data and the parameters are given in Table [H As already mentioned, for a monofractal 
series r(g) is linear with g, and any kind of nontrivial dynamics present in the data 
is reflected in the form of a nonlinear variation of r(g) with g. The signature of 
multifractality in all the original series is clearly visible from their nonlinear variation 
against of r(g) order number g. In some cases the shuffled series as well as the 
surrogated series also show multifractality as strong as the original ones. Note that 
the observable r(g) is directly calculated from h{q), therefore we do not put any 
additional emphasis on it. However, we would like to mention that the GBM model 
reproduces the empirical r(g) spectra very well. 

The most important results of a multifractal analysis is the singularity spectrum. 
To get a quantitative idea about the strength of multifractality, we calculate singu¬ 
larity spectra f{a) for the analyzed series and plot them in Fig. [7] against the Holder 
exponent a. From the top to the bottom panel the diagrams represent, respectively 
the CPI time series, the price time series in China, India and Turkey. The solid lines 
in this figure also indicate the GBM series with the parameter values as given Table 
[T]). As we know that the width of singularity spectrum is a measure of multifractality, 
the multifractal nature of all the original series is well reflected by the corresponding 
singularity spectrum. All the shuffled series within the time period 1993-July 2013 
show much weaker multifractality than the corresponding original series. The same is 
true for the period-I series as well. But for the period-H shuffled series the singularity 
spectral widths are almost as equal as to their original counterparts. The observa¬ 
tions here are consistent with our autocorrelation analysis that is, the time series for 
the period 2010-July 2013 does not capture long-range correlation, and hence the ob- 
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served multifractal nature of these series arise purely from the probability distribution 
functions of their indices. The shuffled series generated spectra, as it is expected, are 
all centered around a ~ 0.5. In some cases (especially Fig. [71(j) and (k)) the peak 
position of the /(a) spectrum is obtained at a slightly higher value of a. A similar 
observation is also made in [52] , where it is shown that the peak position as well as the 
width of a multifractal spectrum may vary with its length along with its type. The 
surrogated series, irrespective of its time period, reproduce stable multifractal spectra 
but much narrower than the corresponding empirical ones. The strength parameter 
of multifractality (Aa) is measured as the width of singularity spectrum at /(a) = 0. 
Analytically that corresponds to Aa = (Ina — ln&)/ln2. Using the (a, b) parameter 
values we calculate the Aa values and are presented in the extreme right column 
of Table [TJ Here we will get a clear and quantitative difference among the various 
series studied. As can be seen in the table that the strength of mult ifr act ality for the 
original long and period-I series is always about two times that of the period-II series, 
while the difference between Aas of the original and the shuffled/surrogated series 
(whenever possible to fit the GBM model) for the series of period-II is marginal. The 
observations indicate that, except the market series in China and India for period-II, 
the observed multifractality in the analyzed gold price indices are partly due to the 
long-range time correlation and partly due to the fat-tailed distribution function of 
the values, and in the exceptional cases it is mainly due to the probability density 
function of the indices; for which we expect Aa (original) ~ Aa (shuffled). The sur¬ 
rogated data generated /(a) spectra in all cases but the series of Turkey are shifted 
towards the lower a side than their shuffled one, but they produce Aa value almost 
identical to that of the shuffled series. In case of the market of Turkey we find no 
change in the day-to-day market return index in several places, and some times the 
index remains constant over a wide section of the series. That might be a reason for 
the exceptional behavior of the market of Turkey. 

4 Conclusion 

Multifractal nature of (i) the consumer price index time series of gold and (ii) the gold 
market time series of China, Indian and Turkey during the period 1993-JuIy 2013, 
has been investigated in terms of the multifractal detrended fluctuation analysis. For 
better understanding of the rapidly increasing and fluctuating market trend over the 
last few years we divide each series into two: first section from 1993 to 2010 (period-I) 
and second section from 2010 to July 2013 (period-II). That is, all together 12 original 
series, corresponding to each of them a randomly shuffled and a surrogated series, are 
analyzed. Multifractal observables, such as the generalized Hurst exponents {h{q)), 
multifractal exponent (r(q')) and singularity spectrum (/(a)) for all the series, are 
extracted and are fitted (whenever possible) to the generalized binomial multifractal 
model series. The following conclusions can be drawn from the our analysis. 

The MF-DFA fluctuation functions for all the analyzed time series nicely follow 
the scaling law ([6|), as is expected for a multifractal series. The generalized Hurst 
exponent spectra corresponding to the original series are found to be order dependent 
and are fitted more or less by the CBM model. The shuffled series as well as the AAFT 
surrogated series also produce order dependent h{q) spectra, but the order dependence 
is much weaker than the original series. The nature of the h{q) spectra are consistent 
with the fact that the multifractality in the long-term CPI time series as well as in 
the market trend of India is generated from both the long-range temporal correlation 
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and the probability distribution of the series values. In the case of the market series of 
China and Turkey the origin of multifractality is found to be the fat-tailed probability 
distribution function, and hence the binomial series f[T5|l cannot describe those series. 
The observed multifractality for period-II is mainly due to the probability density 
function. The r(g) exponent diagrams reflect identical features of the data as the 
h{q) diagrams do. The multifractal singularity spectra for all the original long series 
are significantly wider than the corresponding shuffled and/or surrogated spectra. For 
period-II the original and shuffled series produce almost identical singularity spectra 
that strengthen the observations of the h{q) spectra. 

It is obvious that the market series and the CPI series are highly correlated. 
Therefore, a cross-correlation study [52] between the individual market series and 
market-CPI series may reveal more insight into the gold market pattern. 
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Figure 4: The MF-DFA2 fluctuation functions Fq{s) for the consumer price index 
series of gold plotted with scale s for three different time duration: (a) 1993-July 
2013, (b) 1993-2010 and (c) 2010-July 2013. Predictions from the shuffled series 
(middle panel) and the surrogated series (bottom panel) corresponding to the original 
ones are also shown. 
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Figure 5: The generalized Hurst exponent h{q) spectra for the analyzed time series. 
The lines correspond to the GBM series fllSp . yielding Eqn. fll6p . The best £t param¬ 
eter values are quoted in Table [H 
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Figure 6: The multifractal r(g) exponent spectra for the analysed time series. The 
lines correspond to the GBM series flTHl) . yielding Eqn. (Ha with the same set of 
parameters as optimized for the h{q) versus q data points (see Tabled]). 
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Figure 7: The singularity spectra for the analyzed time series. The lines correspond 
to the GBM series prediction with the parameter values optimized for the h{q) versus 
q data points. The multifractal strength parameters values are given in Table [U 
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